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Abstract 

Pattern formation of sound is predicted in a driven resonator where subharmonic generation takes 
place. A model allowing for diffraction of the fields (large-aspect ratio limit) is derived by means 
of the multiple scale expansions technique. An analysis of the solutions and its stability against 
space-dependent perturbations is performed in detail considering the distinctive peculiarities of 
the acoustical system. Numerical integration confirm the analytical predictions, and shows the 
possibility of patterns in the form of stripes and squares. 
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I. INTRODUCTION 



The topic of pattern formation, or the spontaneous emergence of ordered structures, 
is nowadays an active field of research in many areas of nonlinear science [lj. Pattern 
formation is commonly observed in large aspect ratio nonlinear systems which are driven far 
from the equilibrium state by an external input. Transverse modes may become unstable 
when the amplitude of the external input reaches a critical threshold value, large enough 
to overcome the losses produced by dissipative processes in the system, and a symmetry 
breaking transition developes, carrying the system from an initially homogeneous to an 
inhomogeneous state, usually with spatial periodicity. 

Parametrically driven systems offer many examples of spontaneous pattern formation. 
For example, parametric excitation of surface waves by a vertical shake (Faraday instability) 
in fluids |2j and granular layers [3||, spin waves in ferrites and ferromagnets and Langmuir 
waves in plasmas parametrically driven by a microwave field or the optical parametric 
oscillator 0, 0] have been studied. 

In nonlinear acoustics, a phenomenon which belongs to the class of the previous examples 
is the parametric sound amplification. It consists in the resonant interaction of a triad of 
sound waves with frequencies uj$ , uj\ and LO2 , for which the following energy and momentum 
conservation conditions are fulfilled: 

k = k 1 + k 2 + Ak, (1) 

where Ak is a small phase mismatch. The process is initiated by an input pumping wave of 
frequency uq which, due to the coupling to the nonlinear medium, generates a pair of waves 
with frequencies uj\ and U2. When the wave interaction occurs in a resonator, a thresh- 
old value for the input amplitude is required, and the process is called parametric sound 
generation. In acoustics, this process has been described before by several authors under 
different conditions, either theoretical and experimentally. In the one dimensional 

case (colinearly propagating waves) is considered. In the problem of interaction between 
concrete resonator modes, with a given transverse structure, is studied. In both cases, small 
aspect ratio resonators containing liquid and gas respectively are considered. More recently, 
parametric interaction in a large aspect ratio resonator filled with superfluid He 4 has been 
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investigated [ill ]. 

It is well known that optical and acoustical waves share many common phenomena, an 
analogy which sometimes can be extended to the nonlinear regime [12]. In particular, the 
phenomenon of parametric sound generation is analogous to the optical parametric oscil- 
lation in nonlinear optics. However, an important difference between acoustics and optics 
is the absence of dispersion in the former. In a nondispersive medium, all the harmonics 
of each initial monochromatic wave propagate synchronously. As a consequence, the spec- 
trum broadens during propagation and the energy is continuously pumped into the higher 
harmonics, leading to waveform distortion and eventually to the formation of shock fronts. 
On the contrary, dispersion allows that only few modes, those satisfying given synchronism 
conditions, participate effectively in the interaction process. 

In acoustics, the presence of higher harmonics can be avoided by different means. One 
method is based in the introduction of some dispersion mechanism. In finite geometries, such 
as waveguides [3] or resonators jl^ , the dispersion is introduced by the lateral boundaries. 
Different resonance modes, propagating at different angles, propagate with different effective 
phase velocities. Other proposed methods are, for example, the inclusion of media with 
selective absorption, in which selected spectral components experience strong losses and may 



be removed from the wave field 
dependent complex impedance 



] , or resonators where the end walls present a frequency- 
In this case, the resonance modes of the resonator are 
not integrally related, and by proper adjustment of the resonator parameters one can get 
that only few modes, those lying close enough to a cavity resonance, reach a signifficant 
amplitude. In any of these spectral approach to the problem, in terms of few 

interacting modes, is justified. 

Therefore the selective effect of the resonator allows to reduce the study of parametric 
sound generation to the interaction of three field modes, correspondig to the driving (fun- 
damental) and subharmonic frequencies, and to describe this interaction through a small 
set of nonlinear coupled differential equations. In the present work we concentrate on the 
particular degenerate case of subharmonic generation, where u>i = u>2 and consequently 
ujq = 2o>i, being ujq the fundamental and u)\ generated subharmonic, both quasi- resonant 
with a corresponding resonance mode. This degenerate case has been considered in previous 
experimental studies jjlQ], in the case of small aspect ratio cavities where transverse spatial 
evolution is absent. 
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The aim of the paper is twofold. On one side, a rigurous derivation of the dynamical 
model describing the parametric interaction of acoustic waves in a large aspect ratio cav- 
ity is presented. The derived model is isomorphous to the system of equations describing 
parametric oscillation in an optical resonator, and consequently their solutions are known. 
On the other side, the model predicts the existence of pattern forming or modulational 
instabilities, when the subharmonic field is negatively detuned with respect to the cavity. 
The second aim of the paper is therefore to determine the conditions under which pattern 
formation of sound could be observed in an acoustic system. Numerical integrations under 
realistic acoustical parameters confirm the predicted results. 



II. DERIVATION OF THE MODEL 



A. Three wave interaction in an acoustic resonator 



The physical system we consider in this paper is an acoustic cavity (resonator) composed 
by two parallell solid walls, with thicknesses D and H separated a distance L, containing a 
fluid medium inside, as described in Fig. 1. The different media are acoustically character- 
ized by its density p and the propagation velocity of the sound wave, c. One of the walls 
vibrates at a frequency close to one of the normal modes of the cavity. 

The _ce modes / (eigeufcequeucies) of such resonate, can be calculated by using 
the equation [9( 

( f f \ ( f f \ f 

TZ I tan — — it + tan — — tt J + I 1 — TZ 2 tan — — tt tan — — it J tan — it = 0, (2) 

V Id Ih J V Id Ih J II 

where TZ =p w c w /pc is the ratio of wall to medium acoustic impedances, fn = c w /2D, fn = 
c w /2H, and fi = c/2L are the fundamental resonance frequencies of each individual region, 
respectively. From the numerical solutions of Eq. (j2J) results a non-equidistant spectrum, 
the position of the different modes being determined by the properties and dimensions of the 
different elements. Note that the particular case corresponding to an equidistant spectrum 
(constant free spectral range) is obtained as a limit case when we impose infinite reflectance 
at the walls (TZ — > oo) with negligible thickness. In this case Eq. (0) reduces to tan kL = 0, 
and the modes obey the Fabry-Perot condition k = utt/L. In such a perfect resonator, any 
harmonic of a resonant driving wave is also resonant with a higher-order cavity mode, and 
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the energy flow into these modes leads to wave distortion and invalidates a modal description 
of the problem in terms of the interaction among few waves. 

In a loosy resonator, however, one can get the second harmonic of the driving wave to 
be more detuned than subharmonics with respect to a cavity resonance, thus reducing the 
effectivity of the cascade process into the higher harmonics. This effect is enhanced in the 
case of viscous media, in which the higher frequencies experience stronger losses (absorption). 
Furthermore, it is possible to get subharmonic generation slightly detuned from a cavity 
resonance, which is a necessary condition for the developement of spatial instabilities, as 
will be discussed in the following sections. These two facts justify the description of high 
intensity acoustic waves in a resonator in terms of the interaction among few frequency 
components. Several experimental results demonstrate this fact HQS , and support the 
validity of this assumption in the theoretical approach. 

The main novelty of this work with respect to previous studies, is to consider the diffrac- 
tion of the waves inside the cavity. Diffraction can play an important role when the cavity 
has a large Fresnel number, defined as F = a 2 /XL, where a is the characteristic transverse 
size of the cavity (for example, a 2 is the area of a plane radiator), A is the wavelength and L 
is the length of the cavity in the direction of propagation, considered the longitudinal axis 
of the cavity. Sometimes the case of large F is called the large aspect ratio limit. All these 
assumptions will be taken into account in the derivation of the model in the next section. 



B. Hydrodynamic equations for sound waves 

As a starting point of the analysis, we consider the basic hydrodinamic equations de- 
scribing the propagation of sound waves in liquids and gases, namely the continuity (mass 
conservation) equation, 

^ + V(pu) = 0, (3) 

and the Euler (momentum conservation) equation, 

p (j^ + uv) u = -Vp + /iV 2 u + (p B + |)V(Vu), (4) 

where p is the density of the medium, u is the fluid particle velocity, p is the thermodinamic 
pressure, and p and pb represent shear and bulk viscosities, respectively. Equations (jHJ) and 
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(HJ) must be complemented by the equation of state p = p(p). If the losses due to viscosity 
are small (due just to heat conduction) the process can be assumed to be adiabatic. Then 
the pressure in the state equation can be expanded around the equilibrium and the equation 
of state takes the form 



*2. ^ 2 



p + clp' + ^Tp' 2 + ..., (5) 



where p' = p — po, being p the equilibrium value of the density, Co = a/ (dp/dp) s is the (low 
amplitude) sound velocity and 

oy\ _ 4B 

dp*)~ Pa A' 

where B/A is commonly used in acoustics as the nonlinearity parameter, and has been 
measured in different media Q|. The subscript s denotes the adiabatic character of the 
process, and the ellipsis in Eq. (J5J) the nonlinearities higher than quadratic, which are 
neglected. 

Substitution of Eq. (pj) in (@J) leads to 

p + uV) u = -c 2 Vp' - irVp' 2 +/iV 2 u + (p B + |)V(Vu). (6) 

which together with Eq. (JHJ) are a two-variable model. It is convenient to write Eqs. (JHJ) 
and © in nondimensional form, adopting the following normalizations: 

where V is a reference velocity, small compared with c . Also, time and space are defined 
as 

t = out, x = fcx. (8) 

where u and k are the angular frequency and wave number of a reference wave, and obey 
uj = kco- With this normalization Eqs. (J3J) and (JHJ) have the form 

^ + Mp Vv + Mv Vp = 0, (9) 



d 1 
Mp-^ + M 2 pvV v = -Vp - -FV (p - l) 2 + p M 2 V 2 v. (10) 
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where the losses p, the nonlinearity T and the acoustic Mach number M, are parameters 
defined as 

r = % r - f , (12) 

M = — . (13) 

Co 

In Eq. ()10|) we have used the identity V(Vv) =V 2 v + V x V x v, where the second 
(vorticity) term has been neglected, since its magnitude decays exponentially away from the 
boundaries 



C. Perturbative expansion in the small Mach number limit 

Under usual conditions, the acoustic Mach number take small values (M < 10~ 3 ), which 
allows to treat Eqs. (jHJ) and (0) by perturbative techniques. Thus we consider a smallness 
parameter e as the Mach number, and express the parameters and variables in terms of it. 

Let us assume that in the dispersive resonator, the changes in the shape of the wave as 
a consequence of dissipation and nonlinearity, both along the direction of propagation and 
transverse to it, are small. Also, we take into account that the changes along the transverse 
direction to the propagation, due to diffraction, take place faster than along this propagation 
direction These assumptions allow to consider the problem in terms of fast and slow 
scales. A choice of scales accounting for these changes is 



t = T + er, (14a) 

(z, x, y) = (z, yfex % y/ey) , (14b) 

and expand the state variables p and v = (v x , v y , v z ) as 

p = 1 + epi + e 2 p 2 , (15a) 

v z = v lz + ev 2z , (15b) 

v x = y/ev lx + Ey/ev 2x , (15c) 

Vy = \fev V y + E\fev 2y . (15d) 
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where the order of the transverse components of the velocity is determined by the (slow) 



divergence of the beam 17[. Note that, since at equilibrium fluid is at rest and, on other 
hand, v is nondimensionalized by V (a reference velocity smaller than cq), the variation of 
v is order 0(1). Substituting these expressions into Eqs. (jHJ) and (JTUJ), we obtain equations 
at different orders which can be recursively solved. The leading order 0(e) reads 

dv lz dpi 



dT ' dz ' 
dpi dv lz 



(16a) 
(16b) 



dT dz ' 

which leads to linear wave equations for the particle velocity and density. In general, the 
frequencies of the waves are not resonant with a cavity eigenmode, and the waves at any 
frequency are detuned, the detuning parameter being defined as 

6 t = u\ - Ui, (17) 

where u>i is the frequency of the field and uf is frequency of the cavity eigenmode closest to 
oji. In this case, the solution of the order 0(e) takes the general form of a superposition of 
standing waves 

2 

viz = )A n (x, y, t) sm[u n T - (0 n (r) - 5,T)] sin (k n z) . (18) 

n=0 

where u n = k n . In Eq. ffT8|) we have considered the frequency-selective effect of the walls 
discussed in the introduction, and assume that only three modes, those with frequencies 
obeying u)\ -\- 002 = ^0, can reach signifficant amplitudes. Also, from Eqs. (fT7)|) we obtain 
that 

2 

Pi = ^A„(x, y, t) cos[^„T - (0„(r) - 5iT)] cos (k n z) . (19) 

n=0 

At order 0(e 3 ^ 2 ), the equations for the fast evolution of the transverse velocity compo- 
nents are obtained, 

; (20a) 

(20b) 



dvi x 


dpi 


dT 


dx 


dvi y _ 


dpi 


dT 


dy 



s 



At order 0(e 2 ) we get 

dv u . dv 2z dv lz ld 2 vf z tdp\ _d 2 v lz dp 2 . , 

Pl ^T + ^T + ^ + 2^ + 2^-^ + ^ = °> (21a) 

— + — + + +v — + + = (21b) 

dT dr dz dz z dz dx dy 

Finally Eqs. (J2TJ), making use of the relations in Eqs. (J20J), can be reduced to a single 
wave equation for the density 

d 2 p 2 d 2 p 2 _ 2 d2 viz - d 3 v lz | T d 2 pl | d 2 v\ z | d 2 pi | d 2 pi 2 
dT 2 dz 2 dzdr ^ dz 3 2 dz 2 dz 2 dx 2 dy 2 

where the smaller order solutions appear at the right-hand side as source terms. 

A closed set of equations for the slowly varying envelopes of mode amplitudes A n and 
phases n can be obtained by substituting Eqs. ()18j) and ()19j) in Eq. (|22|). and imposing the 
absence of secular terms in the resulting equation, i.e. neglecting the source contributions 
that contain the same frequency components as the natural frequency of the left-hand side 
part in Eq. (|22j). Otherwise, second order solutions would grow linearly, violating the 
smallness condition ep 2 « p\ assumed in the perturbation expansion. After some algebra 
the secular terms reduce to the following real system for the amplitudes and phases 

dA- 1 

+ (TSiU}iAjA k sin + -puj 2 Ai = 0, (23a) 

deb 1 
Ai-P-- vuJiAjAk cos <p - —V 2 ± Ai - 5iAi = 0, (23b) 
or 2u>i 

where k) = (0, 1, 2), and the other two equations are obtained by cyclic permutations. 
The sign operator Si = +1 for i — 1,2 and — 1 for i = 0. A global phase is defined as 
= 0i + 02 — 0O) stems for the Laplacian operator acting on the transverse space 
r _i_ = {x, y), and a is a coupling parameter defined as 

-H i+ £>- (24) 

Note that in Eqs. (}2*3*j) the terms proportional to p accounts only for viscous losses. 
There are in fact other loss mechanisms, mainly related with finite reflectance of the walls 
or diffraction losses through the open sides of the resonator. When the losses are sufficiently 
small, one can generalize Eqs. (J23|) and consider an effective (phenomenological) loss pa- 
rameter 7j for each mode, just replacing ^uj 2 A n by ^ n A n . The value of these coefficents 
can be obtained experimentally for a particular resonator by small amplitude measurements 



of the decay rate of a given mode, since under this condition (neglecting nonlinearity) the 
amplitudes obey dA n /dr = —^ n A n . 

Finally, for a dissipative resonator, an external source must be provided in order to 
compensate the losses. Consider that a plane wave of amplitude E and frequency uJq is 
injected in the resonator. In each roundtrip, the amplitude of the standing wave will increase 
by 2E, and then 

AA _ IE _c 

where L is the length of the resonator and Cq/L corresponds to the time taken for a wave 
to travel across it. By changing to the nondimensional notation in slow time scale, and 
assuming small amplitude changes during a roundtrip, one can consider the differential 
limit of Eq. (|25jl and incorporate it to the evolution equation of Aq as a driving term. 

A particular case corresponds to degenerate interaction, where U\ = 002- This process de- 
scribes the subharmonic parametric generation. In this situation, the pump and subharmonic 
waves obey the relation uq = 2u\. It is worth to express the resulting system of equations 
in complex form, by defining the complex amplitudes as B n (r±, r) = A n (r_|_, r) exp [i0 n (r)] . 
At this point we also return to dimensional (physical) variables. Since amplitudes A n corre- 
spond, e.g. to dimensionless densities, and pressure is related to density by Eq. ©, then, in 
the degenerate limit we obtain the following coupled equations for the evolutions of pressure 



■flT = -(7o + ido)Po ~ 1 0P1 + l n — V _J°o + -rPin, (26a) 

at po c o ^0 L> 



-7TT = - 7i + WiJPi - * — 2P1P0 + ijj— Vjj?i. 26b 
at PqCq 2u}\ 



together with their complex conjugate. In Eqs. (J26J), pi corresponds to deviations with 
respect to equilibrium pressure values. 

Equations can be further simplified by adopting the following normalizations 

.2p Co7i^ , Q7 \ 

Po = 1 V , (27a) 

auj 

= gogjV^ (27b) 

m**™ ,, (27c) 
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and introducing the dimensionless detuning parameter A n = 5 n / , ~f n . The final form of the 
model reads 

( ' " -(1 + tA )V - Vl + ta V 2 ± V + £, (28a) 



7o dt 

( " 1 -(1 + iA 1 )V 1 + V*V + ia^lVx. (28b) 



7i <9t 

where a n = cf ) /2LJ n 'y n are the diffraction coefficients. This form of the equations is relevant 
for our purposes, since their solutions and stability have been discussed in the context of 
nonlinear optics, after the model given by Eqs. (|2*Hj) has been derived for the degenerate 
optical parametric oscillator. A detailed analysis of the spatio-temporal dynamics of Eqs. 
(j2Hj) has been carried out during the last decade (starting with the seminal work j^), and 
a recent overview can be found in j^j]. In the following sections, we review the basic 
results regarding their homogeneous solutions and their stability, and we study numerically 
the spatio-temporal dynamics under conditions corresponding to a real acoustical resonator. 
Note that the connection between the generic model ()28|) and the particular acoustic problem 
is given by the normalizations performed in Eqs. ()27|) . 

III. MODULATION AL INSTABILITIES OF HOMOGENEOUS SOLUTIONS 

Two stationary states are solution of Eqs. (J28|) : the simplest, trivial solution, 

P °=(TW^ = ' (29) 

characterized by a null value of the subharmonic field inside the resonator, and the nontrivial 
solution 

\V \ 2 = 1 + Aj, (30a) 



|Px| 2 = -1 + A A 1 ± y/ \S\ 2 - (A + Ax) 2 , (30b) 

in which both the pump and the subharmonic fields have a nonzero amplitude, and exists 
above a given (threshold) pump value £ = £ t h- At this value, given by 



\£ th \ = y/(l + A 2 )(l + Al), (31) 

the trivial solution loose its stability and bifurcates into the nontrivial one. The emergence 
of this finite amplitude solution corresponds to the process of subharmonic generation. Note 
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that the fundamental amplitude above the threshold is independent of the value of the 
injected pump, which means that all the energy is transferred to the subharmonic wave. 



These results have been confirmed experimentally for an acoustical resonator in 10 1. 
The character of the bifurcation depends on the detuning values. As demonstrated in 10], 
and also in the optical context, [l9( the bifurcation is supercritical when AqAx < 1, and 
subcritical when AoAi > 1. In the latter case, both trivial and finite amplitude solutions 
can coexist for given sets of the parameters, which results in a regime of bistability between 
different solutions. 

In order to study the stability of the trivial solution (j2HJ) against space-dependent pertur- 
bations, consider a deviation of this state, given by A,- (rj_, t) = Aj + 5Aj (r±, t) . Assuming 
that the deviations are small with respect to the stationary values, one can substitute the 
perturbed solution in Eqs. (j28|) and linearize the resulting system in the perturbations 5Aj. 
The generic solutions of the linear system are of the form 

(SAj, 5 A*) oc e A(fcj - )t e <k - L - r - L , (32) 

where A(kj_) represents the growth rate of the perturbations, and kj_ is the transverse 
component of the wavevector, which in a two-dimensional geometry obeys the relation 
|kj_| 2 = k\ + ky. The growth rates, which depend on the wavenumber of the perturbations, 
are obtained as the eigenvalues of the linear system. This analysis has been performed before 
Q], and we present here the main conclusions, omitting details. 

The eigenvalue (and consequently the instabilities) presents a different character depend- 
ing on the sign of the subharmonic detuning. If Ai > 0, which corresponds to a subharmonic 
frequency smaller than that of the closest cavity mode, the eigenvalue shows a maximum 
at k± = 0, the emerging solution being homogeneous in transverse space, with amplitude 
given by Eq. f|HTT|) . On the contrary, in the opposite case Ai < 0, which corresponds to field 
frequencies larger than the nearest cavity mode, the maximum of the eigenvalue occurs for 
perturbations with transverse wavenumber 



*x = f (33) 

The emerging solution in this case is of the form Eq. (|3*2*|) . which represents a plane 
wave tilted with respect to the cavity axis. This solution presents spatial variations in the 
transverse plane, and consequently pattern formation is expected to occur. 
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Since k± is the modulus of the wavevector, the linear stability analysis in two dimensions 
predicts that a continuum of modes within a circular annulus (centered on a critical circle 
at |kj_| = kj_ in (k x , k y ) space) grows simultaneously as the pump increases above a critical 
value. This double infinite degeneracy of spatial modes (degenerate along a radial line from 
the origin and orientational degeneracy) allows, in principle, arbitrary structures in two 
dimensions. 

The threshold for pattern formation follows also from the eigenvalue, and is given by 



The predictions of the stability analysis correspond to the linear stage of the evolution, 
where the subharmonic field amplitude is small enough to be considered a perturbation of 
the trivial state. The analytical study of the further evolution would require a nonlinear 
stability analysis, not given here. Instead, in the next section we perform the numerical 
integration of Eqs. (|28j) . where predictions of the acoustic subharmonic field in the linear and 
nonlinear regime are given. 

IV. NUMERICAL RESULTS IN THE ACOUSTICAL CASE 

The analytical predictions of the linear stability analysis have been numerically confirmed 
for Eqs. (|2*8Jl in previous studies, in the context of nonlinear optics. In this section we 
demonstrate the adequacy of these resuls for the acoustical case. For this aim, we first 
evaluate the different parameters appearing in Eqs. (|2fi|) for a concrete case. 

Consider a resonator composed by two identical walls of thickness D = H = 0.5 cm made 
of a lead zirconate titanate (PZT) piezoelectric material (q = 4400 m/s, pt = 7700 kg/m 3 ), 
containing water (c m = 1480 m/s, pt = 1000 kg/m 3 ). For this case 1Z = 22.89. The length 
of the medium L can be varied in order to modify detunings. If the resonator is driven at a 
frequency fo = 4 MHz, then subharmonic generation is expected to occur at f\ = 2 MHz. 
The corresponding detunings have been numerically evaluated from Eq. (J2J). For a cavity 
length L = 3 cm, the pump is almost resonant with a cavity mode, Sfo = f{j — fo ~ KHz, 
and the subharmonic is detuned by 5f\ = f^ — fi ~ —1.6 KHz. Furthermore, under these 
conditions the second harmonic at /2 = 8 MHz is highly detuned, by Sf2 = /| — fi ~ —3.7 
KHz, and therefore it will reach a small amplitude. 




(34) 
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The loss coefficients, as stated before, can be obtained for small amplitude measurements 
of the decay rate of each mode in the resonator. In particular, a measurement of the quality 
factor for the different cavity modes, defined as Qi = uji/2^i where performed in [9j for a 
similar interferometer. For the frequencies of interest, the measured quality factors take 
values of the order of 10 3 — 10 4 . From this results we can conclude that reasonable values for 
the decay rates are 70 = 7i = 5 10 3 rad/s, which allows to evaluate the rest of parameters 
in the model. The normalized detuning parameters corresponding to this case are Ao = 
and Ax = —1.6, and diffraction coefficients result ao = 8.7 10 -6 and ai = 2a^ Finally, the 
nonlinearity parameter of water at 20°C is B/A = 5, which substituted in Eq. (|24jl gives 
a = 0.875. 

The theory of the previous section predicts that, when the threshold value Eq. (|34j) is 
achieved, a periodic pattern with a characteristic scale given by Eq. develops. For 

the above conditions, the normalized threshold value is £ p = 1, and the corresponding input 
pressure at the driving wall is obtained from the last of Eqs. (f27|l . and results pi n « 0.1 MPa. 
Also, the wavelength of the pattern is obtained from Aj_ = 27f/k± ~ 1.5 cm. Then, in order 
to observe the pattern the transverse section must contain several wavelengths, which in 
turn which implies a transverse size of 10 cm or more. All these values can be considered as 
realistic. 

In order to check the analytical predictions, we integrated numerically Eqs.([28|) by using 
the split-step technique on a 

linear (pump, losses and detuning) and nonlinear, are calculated in the space domain, while 
nonlocal terms (diffractions) are evaluated in the spatial wavevector (spectral) domain. A 
Fast Fourier Transform (FFT) is used to shift from spatial to spectral domains in every time 
step. Periodic boundary conditions are used. 

As a initial contition, a noisy spatial distribution is considered, and the parameters are 
those discussed above, for which a pattern forming instability is predicted. Figure 2 shows 
the result of the numerical integration. In Figs. 2(a) and (b) several snapshots of the evolu- 
tion at different times are shown, which eventually result in a final stable one-dimensional 
pattern in the form of stripes, shown in Fig. 2(c). 

Numerical simulations for different detunings have been performed. A systematic study 
shows that in most of the cases the system developes stripped patterns with arbitrary orien- 
tations. However, in some cases a pattern with squared symmetry result as the final stable 
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state. An example of evolution leading to squared patterns is shown in Fig. (3), obtained 
for Ai = -4, A = -8 and £ = 1.2. 

V. CONCLUSIONS 

The pattern formation properties of an acoustical resonator where subharmonic gener- 
ation takes place are discussed from the theoretical point of view. A model allowing for 
diffraction of the fields (large-aspect ratio limit) is derived by means of the multiple scale 
expansions technique. The obtained model, which is isomorphous to that obtained for the 
optical parametric oscillator, is analized in detail considering the distinctive peculiarities of 
the acoustical system. A typical acoustical configuration is considered, and the predictions 
of the theory are confirmed by numerical integration under realistic conditions. Numerics 
show that transverse patterns in the form of one- dimensional stripes are usually obtained 
as the final stable state, although the system can support also patterns with more complex 
structures, such as squares. 

We finally note that there has been several experimental attempts to verify the pattern 
formation scenario predicted by Eqs. (|28|) in a nonlinear optical resonator Although 
some transverse patterns with different latticelike structures have been observed, the use of 
other cavities different than planar (confocal and concentric) was required. Owing to the 
complete analogy between the model derived in this paper, and the model describing optical 
parametric oscillation in planar resonators, we believe that the acoustical resonator could 
be a good candidate for the experimental observation of transverse patterns in a planar 
(Fabry-Perot) cavity. Experimental work in this direction is in progress. 
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Figure Captions 



Figure 1. Scheme of a three-element acoustic resonator. Each section is acoustically 
characterized by its density an propagation velocity of sound. 

Figure 2. Developement of stripped patterns for A x = —1.6, A = 0, 7 = 71 = 5 10 3 
and £ = 2, as obtained by numerical integration of Eqs. (|28|). The distributions correspond 
to evolution times t = 0.01 s (a) t — 0.1 s (b) and t — 1 s (c). 

Figure 3. Developement of squared patterns for A 1 = —4, A = —8, 70 = 7i = 5 10 3 and 
£ = 1.2, as obtained by numerical integration of Eqs. (f2~%j) . The distributions correspond to 
evolution times t = 0.01 s (a) t = 0.1 s (b) and t = 1 s (c). 
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